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Abstract 

The quantum search algorithm of Chen and Diao, which finds with certainty a 
single target item in an unsorted database, is modified so as to be capable of searching 
for an arbitrary specified number of target items. If the number of targets, vq, is 
a power of four, the new algorithm will with certainty find one of the targets in 
a database of M items using (1/2) (3(iV/^ ) log4 3 - l) » (1/2) (3(iV/i/o) ' 7925 - l) 
oracle calls, where N is the smallest power of four greater than or equal to J\f. If vq 
is not a power of four, the algorithm will, with a probability of at least one-half, find 
one of the targets using no more than (1/2) (9(JV/ 

j,)i°g4 3 _ ij calls, where v is the 

smallest power of four greater than or equal to vq. 



1 Introduction 



Recently Chen and Diao |l| presented a quantum algorithm for searching an unsorted 
database capable of finding, with certainty, a single target item in an A/"-item database after 
2 |~log 4 AT\ iterations of certain unitary operations. ( \x\ denotes the smallest integer greater 
than or equal to x.) Grassl f|| and Tu and Long [|J have given a recursive implementation 
of these unitary operations, and have pointed out that, with this implementation, the 
number of oracle calls required for the j th iteration increases exponentially with j. 

In this paper I present a modification of the algorithm of [|I] for searching an unsorted 
database of M items for vq > 1 target items, provided that the number of targets is 
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known in advance. In Section 2 below I discuss the case of vq equal to a power of four; in 
this case the algorithm will find one of the target items with unit probability. In Section 
3 I discuss the case of vq not equal to a power of four; in this case the algorithm will 
find one of the target items with probability of at least one-half. The number of oracle 
calls required using the recursive implementation is given in Section 4. The notation and 
terminology follow, in general, those of Jl| and 0. 



2 Number of Targets a Power of Four 

Denote the Af items in the database V by w iy i = 1, . . . , Af. Of these items, a total of v Q 
are members of the subset T of target items. An oracle function f(wi) indicates whether 
a selected item is or is not a target: 

f(wi) =1, Wi e T, 

= 0, otherwise. 

If Af is not already a power of four, we embed the database D in a larger database D 
containing additional non-target items such that the total number of items in D is the 
smallest power of four larger than Af: 

D = VU{w M+1 ,...,w N }, (2) 

where 

N = 2 2n , (3) 

n an integer, i.e., 

n=\\og 4 Af], (4) 



so 



N > Af > A/4. (5) 



The above enlargement of the database is as in Jl]. Here, in addition, we embed D in 
a database D which is four times larger still: 

D = D\j{w N+1 ,w N+2 ,...,Wft}, (6) 

where 

N = AN = 2 2 ". (7) 

That is, 

h = n + l. (8) 

All of the additional items not in T> are by definition non-targets, so equation (0) still holds 
and the cardinality of T is still vq. 

For the database to be searched by a quantum computer Q], t ne N items in D are set 
in one-to-one correspondence with the iV computational-basis states \a\a 2 . . . a2n) : 

Wi ^ \a 1 (i)a 2 (i) . . .a 2 n(i)), i = l,...,N (9) 
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where each of the eigenvalues aj(i) is either or 1. The 2n-component vector of a/s 
associated with Wi is termed the symbol of wf 

S(wi) = a 1 {i)a 2 {i) ■ . .a 2ft (i). (10) 

We also define auxiliary symbol functions 

Sj(wi) = a 1 (i)...a j (i), j = 1, . . . ,2ra, 

S 2 n(Wi) = S(Wi). 

It should be emphasized that the correspondence (|S|) is not chosen to make the symbol 
S(wi) a binary representation of the item index i. On the contrary, it is essential for 
what follows that none of the N items in the set D be represented by states such that 
S2(wi) = 00. That is, we require that 

w t eD=> S 2 (wi) ^ 00. (12) 

(We could, for example, establish the correspondence (|S]) so that Wi G D =>- S 2 (Wi) = 11.) 
Condition (|12"D implies 

Wi ET^ S 2 (w l ) ± 00. (13) 

Extending the technique employed in |jj to the case of multiple targets, we select vq of 
the items with auxiliary symbols S 2 (wi) = 00 to be "ground state items." Specifically, the 
u elements of the set G of ground state items, 



G = {w Gi ,w G2 ,...,w Gvq }, (14) 



are those with the symbols 



S(w Gl ) 


= 00. 


. . 000000 


S(w G2 ) 


= 00. 


. . 000001 


S(w Ga ) 


= 00. 


..000010 


S(w G4 ) 


= 00. 


..000011 



(15) 



The rightmost 2p entries in S(w G ) are all l's and constitute a binary representation of 
u Q — 1, where 

2 2p = u . (16) 
We can now define the auxiliary functions 

fjiwi) = 1 if S 2j = 00...00 but Wi ^G, (u) 
= otherwise; j — 1, . . . , h — p, 

and, in terms of these, the auxiliary oracle functions 

F J {w i ) = f{w i )Vf J {w i ). (18) 

(The symbol "V" denotes logical OR.) Note that 

Ff^{ Wi )=f(wi). (19) 
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The starting state for the iteration is the equally-weighted superposition of compu- 
tational basis states obtained from the state |wgi) — |00 . . . 00) by a Walsh- Hadamard 
transformation, 



A 



\ s o) 



(20) 



Ni=i 



Starting from |s ), a total of to? iterations are performed of the transformation 

\s j+1 ) = -Is.ljlSj), j = 0, . . . ,n? - 1, (21) 

where 



to j = n — p. 



The unitary operator Tj in (^) is defined as 



lj = 1 — 2 \ w i)( w i 

i\F j+ i(wi)=l 



(22) 



(23) 



where I is the identity operator. In terms of its action on computational-basis states, 

X j \w i ) = {-l) F ^^\w i ). (24) 

The unitary operator X Sj in ([2"T|) is defined as 

l s .=l-2\ Sj ){ Sj \. (25) 

The proof that, after to? iterations, the resulting state \s n o) is an equally- weighted 
superposition of the vq states Wi G T proceeds by induction. Using (P0|), (|2T|), (|24] ) and 
(HP, we find, for j = 0, 



A? 



E(-i 



A? 



Wi) 



E(-i 



ko) 



1=1 



A/" V<=i 



(26) 



To evaluate the second sum in (PB"D, divide the set of N states into two groups, those for 
which 5*2 (wi) = 00 and those for which Sa(wi) 7^ 00. The first group contains 2 2 ( n_1 ) states, 
of which the 2 2 ( n ~ 1 ) — v states not in G have F 1 (wi) = 1, and the remaining u states in 
G have Fi(^) = (see eqs. (0), (0)). Of the 3 • 2 2 ("~ 1 > states with S 2 ( Wi ) ^ 00, u 
of these have F\{wi) = 1 by virtue of being target states (f(wi) = 1), and the remaining 
3 . 2 2 («-!) - uq have F^wA = 0. So, 



A 



N 



i2n-l 



i=l 



and p6|) reduces to 



i|-Fl(«Ji)=l 



(27) 



(28) 
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We now assume that for some j, 



\ 8j ) = r*+i £ K)> 

i\Fj(wi)=l 



(29) 



and derive the form of \s j+1 ). From (gl]), (||) and (||) 



)-n+j 



i|F j («) 4 )=l \j|Fj(«)i)=l 



• (30) 



The second sum in (|3~0"D can again be evaluated by counting. The items Wi for which 
Fj(wi) = 1 fall into two disjoint groups, those for which fj(wi) = 1, and the elements of T. 
Of the former group, 2 2{jl ~i~ X) — u have F j+ i(wi) = 1 (those with S 2 j+2{wi) = 00 . . . 00— 
recall that the elements of G are not members of {wi\Fk(wi) = 1} for any k), and the 
remaining ^■2 2<<n ~^~ 1 ^ have Fj + \{wi) = 0. As for the elements of T, all uq have Fj + i(wi) = 1. 
Therefore, 



= 2 2 («-i)- 1 , j = -p- 1. 



i|i^-(wi)=l 



Using ([31] ) in (pQj) , we obtain 



(31) 



(32) 



i|Fj+i(iUi)=l 



After applying n° iterations (pT|) to the starting state (pO]), we therefore obtain (keeping 
in mind that F n o{wi) = f{wi)) 



(33) 



A measurement of |s n o) in the computational basis will with certainty yield one of the 
states corresponding to a target item. 



3 Number of Targets Not a Power of Four 

Only a small number of changes are required in the analysis presented above to produce 
an algorithm which will yield one of the target states with a probability greater than 
one-quarter when the number of targets is not a power of four, and which reduces to the 
algorithm of Section 2 when the number of targets is a power of four. All of the definitions 
through the selection of the ground-state items, eq. (|15|), remain applicable. However, the 
integer p defined in ([16]) must be everywhere replaced with p 

2 2p = v, (34) 

where v is the smallest power of four larger than Uq. I.e., 

p = [log 4 z/ ] , (35) 
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v > v > v/ 4. (36) 

The rightmost 2p entries in S(wg^ ) constitute a binary representation of u — 1, but they 
will are not all l's. The definitions (|17]) and ( |l8l) of the auxiliary functions fj{wj) and the 
auxiliary oracle functions Fj(vJi) remain unchanged. However, most significantly, eq. ([19]) 
is replaced with 

{wi\Ffi-p{wi) = 1} D T, (37) 

since not all items with SW-p) = 00 ... 00 are in G. 

So, a derivation parallel to that in Section 2 leads to the conclusion that, by beginning 



with the initial state (EOT) and performing n — p iterations (EHI), we obtain the state 



\Sn-p) = 2 P \ W i)- 
i\F a _p(wi)=l 



(3J 



If a measurement in the computational basis is made of the state (|38D, the probability that 
one of the target states will be obtained is 



where 



Po(p) = P, 



"0 

p = — ■ 

V 



(39) 



(40) 



The probability of finding a target state is thus between one, when u = v (p = 1), and 
somewhat above one-quarter, when i>o = v/4 + 1 (p = l/4+l/V). 

Now suppose that, rather than making a measurement after n—p iterations, we perform 
an "extra" iteration, i.e., compute 



'n— p- 



(41) 



before measuring. The definitions fllTD , (p~8|) of fjiwi) and Fjiwj) work for j > n — p and, 
with the relations ([34|) , (|36l), imply that, regardless of the value of u , 



Fn-p +q (wi) = f(wi), q>l. 



For j = n — p the summation formula corresponding to (|31~D is 

yip , \ 1 f-D^-f+iw 



^i|F fi _p(«)i) 

The state resulting after one extra iteration is 



Ei|F fl _ j3 K)=i(~ 1 ) /(M ° 

2 2 p - 2t/ 



(42) 



(43) 



l«- - A — 9-P+ 1 

l^n— p+1/ ^ 



1-5) |Wi>-5 Yl \ W i) 

i\fi w i) = 1 *|/n-p=l 



where 



5=(4p-l)/2. 

The probability of obtaining a target state upon measuring \sn-p+i) is 

P 1 (p)=p(3-4p) 2 . 



(44) 

(45) 
(46) 
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For 1/4 < p < 1/2, Pi(p) > P (p), while, for 1/2 < p < 1, P x (p) < P (p). So, the 
appropriate strategy is to make a measurement after 



rij=n — p 

iterations ifl/2<p<l, and to make a measurement after 

nf = h — p + 1 



(47) 



(4£ 



iterations if 1/4 < p < 1/2. The probability of obtaining a target state will in this way be 
at least as large as P (l/2) = Pi(l/2) = 1/2 (see Fig. 1). 
Yet another iteration before measurement gives 



| s h-p+2 



> 



-p+1 



where 



i|/(to<)=l i|/ fl _p(u>i)=l 



[l-5) 2 p-5 2 (l-p) 



C = 8 

and a probability of target-finding of 

P 2 (p)=Ap{l-5f(l-C) 2 



(49) 



(50) 



(51) 



Despite the extra iteration, the probability of obtaining a target state when p = 1/2 is not 
increased; P2(l/2) = 1/2. This is true for an arbitrary number of additional iterations. 
The quantum state obtained after n — p + q iterations, q > 1, is of the form 



\ s h-p+q) 



)-P+l 



(52) 



where A g and B q satisfy the recursion relations 



A. 
B, 



9+1 



Al P -Bl{l-p) 



A 



9+1 



1 + 



A 2 n-B 2 (l-p) 



The probability of finding a target upon measurement is 

P q {p) = AA 2 q p. 



(53) 

B 9 . (54) 

(55) 

From ( PI ) and ( [52] ) we see that Ai = 1/2 and B\ = —1/2 when p = 1/2. The relations 
(H)-® then show that 

P q {l/2) = 1/2 V q > 1. (56) 

This is not in any sense to claim that iteration algorithms different than those considered 
here might not improve on the probability of finding a target when p = 1/2. Nor is it to 
say that iterations beyond h — p + 1 necessarily have no use. Probability functions P q (p), 
q > 2, can, for values of p ^ 1/2, be larger than either -Po(p) or -Pi(p), indeed as large as 1 
(see Fig. 1). 
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4 Required Number of Oracle Calls 

Grassl and Tu and Long |J have presented the following implementations of the oper- 
ators Ij and I s ., and have evaluated the number of oracle calls required each time these 
operators are applied. From eq. ( j2"4"|) we see that Ij can be written as 

^ = £(-i) Fi+lK) K><< (57) 

i 

From the condition ( |i~2"|) on the representation of elements of .D (and, therefore, on all 
elements of the target set T), and the definitions (|T7|), ( |T%| ) of fj, Fj, it follows that 

Therefore 

ii= m-2 £ H-^^WKI l ( 59 ) 

V iKeT J \ k J 

and we see that each application of Ij requires a single call to the oracle, since the f/s are 
independent of /. 

From the iteration condition ([21]), the definition ([25]) of J s ., and the unitarity of ij and 
i Sj , we see that the operators I s satisfy the relation 

/ s , +1 (60) 

Let i(j) denote the number of oracle calls required by i Sj . Since ij requires one oracle call, 
( |60|) implies 

t(j + 1) = 3t(j) + 2. (61) 

For j = 0, 

which is independent of /, so 
and t(j) has the closed form 

t(j) = V - 1. (64) 
Taking into account the single oracle call required by Ij, the total number of oracle 



7 ao = /-2|s >(so|, (62) 
t(0) = (63) 



calls required for rij iterations of ([21] ) is 



nj — l 

CM = E + ^ (65) 

i=o 

which, using (|64]), has the value 

C(m) = (1/2) (3^ - 1) . (66) 

It follows from the results of Section 2 that, for Vq a power of four, the required number 
of oracle calls to obtain a target with unit probability is 

C° = CM = (1/2) (3(iV/z/ ) log 4 3 - l) . (67) 
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If vq is not a power of four, the results of Section 3 imply that the number of oracle calls 
to obtain a target state with probability of at least one half is 

C> = C(n>) = (1/2) (3(N/v) l °^ 3 - l) (68) 

if p — Vq/v is between 1/2 and 1, and 

C< = C(nf) = (1/2) (9(N/v) l °^ 3 - l) (69) 

if p is between 1/4 and 1/2. 

The original algorithm of Chen and Diao [|IJ] performs two series of n iterations of ([21]) , 
so the number of oracle calls required to find the unique target item by that method is 

C CD = 2C{n) = 3N l °^ 3 - 1 . (70) 

The exponent log 4 3 is approximately equal to 0.7925. So, with this particular imple- 
mentation of the operators Ij and I Sj , the computational complexity of the algorithms of Jl| 
and the present paper scales more slowly than that of the best possible classical algorithm 
(0(J\f)), but not as slowly as that of Grover's algorithm |5| (0(VJ7)). Unlike Grover's 
algorithm, these algorithms will find a target item with certainty^ if the number of targets 
is a power of four. It is not known at present whether the implementation employed here 
is the most efficient possible, or if implementations requiring fewer oracle calls may exist. 
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Figure Caption 

Figure 1. Probability P q of finding a target with q "extra" iterations, as a function of p. 
Solid line: q — 0. Dashed line: q — 1. Dotted line: q — 2. 
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